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■ld 1 • 

position of pole of orbit in celestial sphere . Long. 352*8 
1 °?! .. Lat. + 64*6 

:§■ - ■ k.a. 316*1 

Deel. + 53*8 
1 


■ 1 
.Si 


hThese elements give for the Mass of Mar.* 


3,090,000 


The Liner Satellite. 

Major semi-axis of apparent orbit at distance [9*5930] . 33"*o± 1" 

Period of revolution ........ 7 h 38 m *5 ± o m *5 

Hourly motion in Areocentric longitude .... 47°*n 

Passage through the eastern apsis (p => 70°) Aug. 20, IS^’O, W.M.T. 

Yery respectfully, your obedient servant, 

Jom* Rodgers, 

Rear Admiral and Superintendent. 
The Hon, R. W. Thompson, Secretary of the LTavy. 


On the Solution of Kepler's Problem . 

By J. W. L. Glaisher, M.A., F.R.S. 

1. Writing Kepler’s equation in the form 


x = a + e sm cc, 


(0 


the solution in powers of e is 


a? = 


= a+e sin a + ££ (sm> a) + j~ ^ (sin* a) + & 0 „ 


and, on performing the differentiations, it will be found that 


x == a 

+ e . sin a 
+ e 2 . sin a cos a 

4- e 3 . sin a cos 2 - sin 8 a 
2 


4- e 4 . sin a cos 3 a — - sin 3 a cos a 
3 

11 . 


13 

+ e 5 .sin a cos 4 —- sin 3 a cos 2 «+ ~ sin 5 a 
3 24 

+ e 8 . sin a cos 5 a—— sin 3 a cos 8 <24- — sin 5 a cos « 

3 15 

4- e\ sin a cos 6 sin 3 $ cos 1 a 4 -sin 5 a cos 2 a— rA sin 7 a, 
6 120 720 


54 i . 7 


+ e 8 ,sin a COB 7 a — — sin 3 a cos 5 <H* ■ - sin- 1 a cos 3 a — ~~r S4n7 a cos a 

3 3 4 o 

+ &c. ....... 


(2) 


K K 2 
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IThus, as is well known, the expression 


from the root of the equation (i) only by quantities of 
the order e 3 . It is also to be noticed that the second term of 
(3) represents the sum of the terms contained in the first column 
of (2), continued to infinity. 

It is evident, therefore, that if the value of the series (2) 
had to be calculated up to the terms of the fourth or a higher order, 
it would be preferable to calculate the value of the second term 
in (3) rather than the values of the terms in the first column 
of (2), so that, in the calculation of the terms in the other 
columns, we may assume that the value of log (1— e cos a ) has 
been already obtained ; and it is of some interest to examine 
whether these terms can be arranged in such a manner that the 
knowledge of the value of log (i — e cos a) may be applied to 
reduce the number of logarithmic entries required for their 
calculation. 

§ 2. It is clear that the terms in e 3 and e 4 in the second 
column of (2) can be represented by 

I e 3 sin 3 a 


» , 2 (1 -e cos a) 


so that 


a + --- -- 

i—e cos cl 2 cos a ) 


e sin a 1 e 3 sin 3 a 


differs from the root of (1) only by quantities of the order e 5 : 
in fact we find that 


e sin a 1 e s sin 3 a 


x « a + 


i—e cos a 2 cos a y$ 



sin 3 a cos 2 a + 




+ c 7 . - 



sm 3 a cos 4 a + 


1291 

120 




+ &c. 


( 5 ) 
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§ 3. The terms involving e 5 and e 6 in (5) 


= e 5 . —- sin 3 a 

18 


+ — sin 5 a 
72 


20 . ^ 

+ e 6 .—rr- sin 3 a cos a + sin 5 a cos a , 

81 


1369 

4°5 


5= c 5 . — * 0 (1 + — e cos a) sin 3 a 


43 . 10952 . . _ 

+ e°. — (1 + —— e cos a) sm 5 a; 

72 1935 

SO that these terms may be represented by 


e 5 sin 3 a 43 
- ; + ' 


sin 5 a 


.18 (\—e cos a) 4 * 0 7 2 (i — e cos 


. . ( 6 ) 


If the exponent IQ 95 2 be replaced by —, the first two terms of 
1935 3 

the expanded value are 

43 e 5 sin 5 a (1 + - 3 - e eos 

72 V 1935 

which differs from the quantity to be represented by only 

- 43 - e 6 sin 5 a cos a = — e 6 sin 5 a cos a , nearly. 

3240 249 


Thus the expression 
where 


a + A, 


A* 


e sm a 


e z sin 3 a 


e 5 sin 3 a ,43 e 5 sin 3 a ^ 

40 *i" , vlj 17/ 


I —6 cos cl 2 (1 __ e cos a )“ 18 (1 —e cos «) 4 9° 7 2 (1 — e cos a) 1 *’ 


differs from the root of the equation (1) by 


■ 43 — 6 e gin 5 ^ cos a + terms of the seventh and higher orders, 

3240 

so that practically (7) gives the solution of Kepler’s problem to 
the sixth order inclusive. 

§ 4. Supposing the values of log e, log sin a , and log cos a 
to be given (in addition, of course, to the logarithms of the 
coefficients, which are assumed to be known), the calculation 
from (2) as it stands, up to the sixth order inclusive, would 
require twelve entries of the table of logarithms, while the 
calculation of log (1 — e cos a) requires two entries, and that of 
the terms in A four more, so that by the use of A six entries are 
saved. If log (1 — e cos cl) be calculated and used in obtaining 
the first column of (2), the other terms in (2) require six entries, 
making nine in all, so that three are saved by the use of the last 
three terms of A. The calculation, of the expression (4), which 
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gives the value of the root up to the fourth order inclusive, 
requires four entries, while the calculation from (2) requires six 
entries if the values of the separate terms in the first column are 
obtained, and five entries if they are replaced by 


e sm a 
i —ecos a 


m § 5 - On developing the expression A up to the seventh order, 
it is found to be 


= e .sm a 


+ e 2 . sin a cos a 


4- e s . sin a cos 2 a — — sin 3 a 


+ e\ sin a cos 3 a — — sin 3 a cos a 


+ <? 5 .sin a cos 4 a — ^ sin 3 a cos 2 a — ^ sin 3 a + —■ sin 3 a 


18 


72 


+ e 6 . sin a cos 5 iz- — sin 3 a cos 3 a— ~ sin 3 a cos a + ~ sin 5 a cos a 
ol 51 210 

+ e 7 . sin a cos 6 sin 3 a cos 4 « ——■ sin 3 a cos 2 a + sin 5 a cos 2 a 

2 43 7 2 9 324 

+ &c. 


e .sm a 

+ e 2 . sin a cos a 


+ 6 3 . sin a cos 2 a - 


sin 3 a 


+ e 4 . sin a cos 3 a — -- sin 3 a cos a 


11. 12’ 

+ e 5 .sin a cos 4 a— — sin 3 « cos 2 a + — sin 5 a 

3 24 

+ e 8 .sm a cot 5 a- 4- sm 3 a cos 3 a+ sm 3 a cos a 

3 648 

+ e 7 .sin a cos 6 sin 3 a cos 4 « + — sin 5 a cos 2 a 

729 2916 

+ &e. 

Subtracting this from (2), it appears that the correction to be 
added to A to give the true value of the root of the equation (1) 


is 


13 

e 6 . —-sin 5 a cos a 


+ e 7 . 


3240 

B 7 S-# 

1458 


sin 3 a cos 4 «.+ sin 6 a cos 2 a- — s ; n 7 a 


29160 


720 


+ &c., 
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! which very nearly 


| = e*. -sin 5 a cos a 

249 

I ° 9 

+ e 7 . - sin 3 a cos 4 a + sin 5 a cos 2 a—— sin 7 , a 

292 20 4 

+ &c., .* ( 8 ) 

The smallness of the coefficients of the terms of the seventh 
order in (8) compared with those in (2) is worthy of remark: 
the sum of the first three coefficients in the e 7 line in (2) is 
equal to about 23 and in (8) to less than unity. 

It is noticeable also how nearly the expression 


I e 3 sin 3 0 


2 (1 —e cos 


represents the first three terms of the second column of (2), the 

coefficient of the third term differing from ~ by only The 

quantity to be added to the expression in (4) to give the value 
of the root up to the fifth order inclusive is 


which 


since 


e 5 . 


1 • 2 

— sin 3 a 
18 


+ —- sin 5 a, 
72 


I 3 

= e 5 . — sin 3 a + ™ sin 5 a , nearly,] 


3 _ 43 _ 1 
5 7 2 3 6o ‘ 


§ 6. The numbers within square brackets denoting the Briggian 
logarithms of the coefficients expressed in seconds, the value of 
A is 


r - e sin a 

[5-3144^513 


-[5'0i3395i] ( - 


e 3 sin 3 a 
— e cos af£ 


[4-0591526]- 


t? sin 3 a 

\4_P 

— e cos a) 9 


+ [5* 0 9 0 5 611 ] 


e 5 sin 5 a 
(1 — e cos a )^ 7 


As an example of the application of this formula, let a = ^ 

and e = 50°, this being the example considered on pp. 376, 377 
of this volume. 

The work is then as follows, every figure that has to be written 
down being printed:— 
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log e = 9*3979400 log cos a = 9-8080675 

,, 6 3 = 8*1938200 „ sin a = 9*8842540 

„ e 5 = 6*9897000 } , sin 3 « = 9*6527620 

sin 5 <2 = 9-4212700 

log 6 9'39794°° 

„ cos a 9-8080675 

9-2060075 
e cos tz 0-1606969 
1 — 6 cos a 0*8393031 


log (I ■ 

— 6 COS <2) 9’9239l89 

Ar. Co, 0-0760811 

= log P, say, 



0008453456 

= log Pa 



0-25360368 

= log P s e° = log V 

LO 

2 


0-33813824 

= log PV° 



0*42267280 

= log P 5 ^ 0 



0-43112626 

= log P's 1 = log P 

17 

3 

(+) 

(-) 

(-) 

( + ) 

5*3*4425* 

5' OI 3395 I 

4’°59 I 5 

5-09056 

939794°° 

8-1938200 

6-98970 

6*98970 

98842540 

9-6527620 

9-65276 

9‘42I27 

0-0760811 

0*2536037 

0-3381:4 

O43II3 

.4-6727002 

3-1135808 

103975 

r93266 

// 

// 

n 

// 

47065-23 

1298-91 

10-96 

85-64 

85-64 

10-96 



47 i 5°-87 

1309-87 



1309-87 




4584I-00 = 

12° 44' i"-oo; 



thus A = 12 0 44' i'^oo, and 




x — a + A, = 62° 

44' i"'oo, . . . 


the true value being (see p. 377) 62 

0 43 ' S6"*oo. 



Tlie value of the expression in (4), that is of a + the first two 
terms of A, = 62° 42' 46'^ and the correction for terms of the fifth 
order, viz. : 

1 g 

e 5 . —« sin 3 a + - sin 5 a 

18 5 

= — 5^03-f 3 i// - 88 ; giving x == 62° 43' 13" 
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LO 1 

^1 

1 ! §7. The correction to be applied to (7) in order to inclnde 
Approximately the terms of the seventh order, viz.: 

«' e 7 .— sin 5 a-— sin 7 a 

| 20 10 

r- 

r- 1 

o // '5o~i // 75 = ~ i //# 25. As the error of (9) is about 5",it will 
~Be seen that, the seventh order correction just obtained being 
very small, there is very little advantage in including these 
terms. In general, the classification of the terms according to 
orders of e supposes that e is very small; so small, in fact, that 
a term involving e n+1 may be regarded as distinctly smaller than 

one involving e n : but if e be as large as -, the terms of the orders 

4 

near to e 8 do not differ much from one another in magnitude, in 
consequence of the increased values of the coefficients, which 
somewhat counteract the effect of the additional factor e. This 
is readily seen by considering the terms in e 8 . By expanding 
the terms in A, it will be found that the terms in e 8 , omitting 
the first, are 


■e — 14*90809 sin 3 a cos 5 a— 1*44389 sin 3 a cos 3 a + 28*82887 sin 5 ct cos 3 a 
—16*35198 sin 3 a eos 5 a + 27*38498 sin 5 a cos 3 a ..(10) 

(in which the coefficients are expressed as decimals so that their 
magnitudes may be more evident): the terms in e 8 in (2), 
omitting the first, are 

e: —16*33333 sin 3 a cos 5 a + 28*46667 sin 5 a cos 3 <2—6*21270 sin 7 a cos a 

.... (II) 

and the difference 

= e 8 . 0*01865 sin 3 a cos 5 a+ 1*08169 sin 5 a cos 3 a — 6*21270 sin 7 a cos a : 


in the case of e = this difference = o //, oo + o //, 24—1 //# 95 

4 

= — l"'*]!. 

The close agreement between the coefficients of sin 3 a cos 5 a 
and sin 5 a cos 3 a in (10) and (11) is remarkable; and it will 
be noted that the coefficient of the new term, viz. the term in 
sin 7 a cos a , is, as in the case of the terms of the seventh order, 
much the largest. 

It thus appears that if e be so large that, in order to obtain 
the value of x to the nearest tenth or hundredth of a second, 
terms of a higher order than the sixth have to be included, in 
general no great additional accuracy will be obtained by applying 
to A the corrections for the terms of the seventh and eighth 
orders, while of course if the terms in e 7 may be neglected, the 
corrections are insensible. Whenever, therefore, the series in (2) 
affords a suitable solution of the problem, the expression .a +A 
gives the solution up to the seventh order, and in general but 
little additional accuracy is gained by applying the corrections 
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i^pi 

j :|for the seventh, and higher orders. The advantage of using the- 
^formula a-{-A would be greatest if e were very small, so that 
;mjhe order of e governed the magnitude of each term, and if the 
llyalue of x were required with great exactness ; as, for example, 
gif the accuracy required were such that ten-figure logarithms 
l”had to be employed in the calculation of the first term of A, but 
" that terms above e 7 were insensible. 

The following are the values of the terms in (2), for this case 

of e == a = 50°, up to the sixth order inclusive :•— 

4 - 

£ = 5°° 

+ 395 o 2 " , oo 

+ 63 47"-8s 

+ io2o"-o8 - 724 // ‘40 
+ i63"*92-388"-03 

+ 26"*34— 13 7"* 18 + 28"* 78 

+ 4 / '*23- 4o"*o8 + 26"*76 

= 5o° + 4583o"*27 = 62° 43' 5o"*27; 

the terms in e 7 

= + o"‘68 - io"*47 +14 //, 76 - i"*46 
== 3 "* 5 U 

and the terms in e 8 

= o"*n- 2 ,/ *54+ 6"'28 -i // ‘95 
= i"’ 9 °; , 

so that, when the seventh order terms are included, the value is 
62° 43' 53 // *78; and when the terms of the eighth order are 
also included, 62° 43' 55 //# 68. 

§ 8. As another example, let e = —, a = 50°, and the terms 
of A 

= i6886 / '-23-57"-85-o"'07 + o" , 47, 

so that a + A — 54 0 .40' .28"78, which will be found to satisfy 
the equation to within the hundredth of a second. 

§ 9. The exponents in the denominators of the terms in 

A are 22 , —, 21 . if the last exponent be taken to be 22 the 
9 9 9 9 

sixth order error 

= + & sin 5 o, cos a, - + e 6 sin 5 a cos a , nearly, 

and is thus much larger than when the exponent is , the error 
in this case (see § 3) being 

__ __L_ e 6 g | n 5 0 cog nearly. 
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w 

I'vPl 

; § 10. If the terms in (2) be treated in columns as in §§ 2 

;%nd 3, each column being kept distinct, it will be found that 
[cijLp to the sixth order inclusiye, 



e sm a 

a + - 

1 —g cos a 


L e3 s ^ n3 a J g5 s ^ n3 a C0§2 a 13 gS sin 5 a 

2 (1 — e cos a ) X 3° 18 (1 — e cos af* 2 4 (1 —e cos a) s ®$ 


so that 


where 


x — a + B, 


^ e sin a I e 3 sin 3 a I c 5 sin 3 a cos 2 a 13 e h sin 5 a 

i — e cos a 2 (1 cos a) 1 * *8 (1 cos a,) 4 * 24 (i cos a) 3 * 3 * 

the correction to be added to B being 


e 6 . sin 5 a cos a 
480 

+ e 7 .—sin 3 a cos 4 a + sin 5 a cos 8 a—— sin 7 a 
1458 3840 720 


+ &c., 

which, very nearly, 


= e 6 . — sin 5 a cos a 

53 


+ e 7 .-^- sin 3 a cos 4 a + — sin 5 a — sin 7 a 
292 4 

+ &c. 


It is clear that the expression B is not so convenient for calcu¬ 
lation as A ; and also the residual sixth order term is larger for 
B than for A. 

§ 11. In t. v. (1859) the Ami ales cle V Observatoire de Paris , 
pp. 349, 350,* M. Serret gives an account of a method, indicated 
by M. Leverrier on p. 192 of t. i., of solving Kepler’s problem 
by means of a series of a different form to (2). The method is 
as follows r— 

Let x 0 be an approximate value of and put 


then 


v = a—x 0 + e sin x 0 ; 


x—e sin x — x 0 — e sin x 0 + v. 


This is the equation to be solved, and we can calculate x by 

* “ Note sur l’equation dont depend l’anomalie excentrique, et sur les series 
qui se presentent dans la theorie du monvement elliptique des corps celestes 
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developing it in ascending powers of v by Maclaurin’s theorem ; 
we have 


f \ dx 

(l-eco SX ) Tv 



. N d 2 x . / dx \ 2 

(*-*<*»*)&* c sm * = o, 

, x d 3 x . dx d 2 x / dx \ 3 

(I-, cos *) -^ + 3 e mg - — + e cos * = o, 

&e. &c. 


When v = o, x ■= a? 0 , and 


whence 

a? = + 


<?? 

tfo I 


dv l — e cos x n 


d 2 x 

dv 2 


e sm x 0 
{i — e cos # 0 ) 3 ’ 

e cos x 0 


d s x 3 e 1 sin 2 x 0 

dv 3 “ (i-e cos #o) 5 _ (1—« cos tf 0 ) 4 ’ 
v _l esinxo 9 


l — e cos x 0 2 (1 — e cos x^f 


f i e 2 sin 3 x 0 i e cos x 0 1 3 
U2 (1— e cos x o y 6 (1 — e cos x o yj + c ' 


(12) 


In the case of an orbit in which e is very small, put x Q = a, then 
v = e sin a and, neglecting terms of the fourth order, 


e sm a 1 

x = « +- 

1 — e cos cc 2 


( e sin a \ 3 
i—e cos a) 


(13) 


M. Serret remarks that this formula is very convenient for 
calculation, and that the last term can be immediately derived 
from the previous term. 

§ 12. It will be found that the coefficient of v A in (12) is 

_ 5. e3 sin 3 x q _5_ i 2 sin x Q cos x Q e sin x 0 
8 {i—e cos x 0 ) 7 12 ( i—e cos x 0 ) 6 + (1 — e cos x 0 f ’ 

and therefore, putting x 0 — a , the formula, whose first three 
terms are given by M. Serret’s equation (13), is 


e sin a i e 3 sin 3 a 

x = a + ------- 

l — e cos a 2 {i—e cos a) 3 

^ y e 5 sin 5 a I e 4 sin 3 a cos a 

2(1-6 cos a) 5 6 (1—6 cos a) 4 J 

_ 5. g 7 ,sin 7 « ^ 5 6 6 sin 5 a cos a 1 e 5 sin 5 a 

8 (1 —e cos a ) 1 12 (1 — e cos a ) 6 + 24 (1 — e cos a) 5 

+ &c.; ..... (14 

the second line being derived from the terms in v 2 , and the third 
from those in v$. 
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I 

1^1 

U; In the case of e = a = 50°,"this gives 

4 

x = 5o 0 + 47o65"-23-i 225"*24 
p; + 63"-79— 78^20 

8! - 4 "‘i 5 + 9 7 , i6 + 5 v, 32 

+ &c. 

o / // 

The first line = 52 43 59-99 

the second line = — 14-41 

the third line = + 10-33 


5 2 43 55-9I 

The close approximation afforded by the first line (the error 
is only 3 f/, gg) is very remarkable, but it is of course accidental; 
the term of the fourth order = 78"*20, and those of the fifth 
order amount to 6g"'ii. 

§ 13. The formula (13), which is true to the third order, only 
differs from (4), which is true to the fourth order, in the 

exponent of the denominator, which is — instead of 3. If 

3 

e =r —, a = 50°, the last term of (13) = 56 // *59 ; and the value 
10 

of x from (13) is 54 0 40' 29 // *64, and from (4) is 54 0 40' 28' / ‘38, 
the true value as found in § 8 being 54 0 40' 28"'78. 

As another example, let e = -, a being, as before, 50°, then 

7 

(14) gives 

% = 50° + 24854 // *92-180"-45 
+ 2 / '-62— 6"-o8 

— o"*o5 + o''-22 + o ff, 22 

+ &c.; 

the first line = 56 6 51' i4' v '47, the second — 3 // ‘46, and the third 
°"*39, leading to x =2 56° 51' n'^o ; and the true value obtained 
by substituting this in (1) is 56° 51' ii"*39. 

The value of a + A is 

50° + 24854"-92- i86"-34-o"-47 + 3"*34> 

which = 56° 51' 

The value obtained from (4) is 56° 51' 8"*58. 

§ 14. It will be noticed that in both these examples the 
formula (13) affords a very close approximation to the true 
value; this is due to the fact that the terms in the second 
line of (14) nearly counteract one another for the values of 
e and a that have been taken; for, these terms 


1 e i sin 8 a , . _ I , 

— -- jg sm z a cos a], 

2 (1 -e cos ay 1 3 / 
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i^nd the term in brackets is small for the assumed values of 
and a. As an example in which this is not the case, let 

'M = 7». ® = I 3 °°) fhen (14) gives 

:i: 4 

1 

l 30 ° + 34032''-99 - 463"'25 
+ 12" -61 + 21 "-38 

— o' ; ' 43 - i"- 46 +i"*c> 5 . 


b The first line = 139° 19' 29"*74, the second 33' / *99, and the 
third — g^‘ 84, leading to x-= 139 0 20' 2"'89, 

The value of a 4-A 


= 130° + 34032 // *99 - 44o ,/ -8 i - 2"*59 4 -13"*64 
= 139 0 20' 3"-23. 


The value of x given by (4) is 139 0 19' 52 //, i 8, which only 
differs from the true result (viz. 139° 20' 2 //> 93) by about 11 ". 
It is evident that the effect of the denominators in (7) and (14) 
is to increase the magnitudes of the terms if a be in the first or 
fourth quadrants, and diminish them if a be in the second or 
third quadrants. 

The expression in (4) is only very slightly more complicated 
than (13), and as it represents the root to one order higher, it is 
in general to be preferred. 

§ 15. The series in (14) is clearly the expansion of x in powers 

of -, and it may therefore be obtained directly in the 

1 — e cos a 't 

following manner:— 

Let 

-e 

g = - , 

i —e cos a 

then 

e 

e ss --• 

I + e cos a 

the equation is 


viz. 


€ 

sc = a+ --— 

1 + e cos a 


sin sc, 


sc = a + e (sin a?—a? cos a + a cos a). 


and therefore, by Lagrange’s theorem, 


sc - a + € sin a + 


e 2 


1.2 


(sin a?—a; cos a + a cos a) 2 


+ 


6 s or 

1.2.3 


(sin a ?—x cos a + a cos a) 3 + &c., 


• (IS) 
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lo 1 

* 

x being put equal to a after the differentiations have been per¬ 
formed ; thus 



r ~- 1 


ooi 



x ~ a 
+ e . sin a 

„ 1 . „ 

+ e 3 . — sin 3 a 
2 

, 1 . „ 

+ e 4 . — — sm J a cos a 
o 

+ &c. 





It is evident that there will be no terms in the expansion in¬ 
volving powers of a cos et, for 

sin x—x cos a + a cos a — sin a , when x = a, 


and its differential coefficient = cos x — cos a . When x = a this 
vanishes, and the terms are in consequence greatly simplified: 
in fact, the development of (15) may be very rapidly effected. 

§ 16. A similar method affords the expansion of x in powers 

of--- ; for, let 

1 + e cos a 

e 

7] =-, 

i + e cos a 

then 


viz. 

and therefore 


v • 

x — a + --— sm x, 

i— 7 ] cos a 

x — a+ 7 ] (sin x + x cos a—a cos a), 


x s= a + 7 ] sin a + 


1.2 


~ (sin x + x cos a—a cos a ) 2 


t? d 2 

+ 1.2.3 dx 2 


(sin x + x cos a —a cos a) z + &c., 


x being put equal to a. 

We thus find 

x » a 

■4-17 , sin a 
+ rf. 2 sin a cos a 
+ 77 s . 4 sin a cos 2 a — sin 3 a 


+ n 4 . 8 sin a cos 3 a—— sin 3 a cos a 
o 

+ 7j 5 .16 sin a cos 4 a — — sin 3 a cos 2 a t 
3 

+ &C. 


— sin’ a 
24 



© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjournals.org/ by guest on March 21, 2015 






Mr. Glaisher , On Solution of Kepler's Problem, xxxvn. 9* 
The first column evidently 


i^ri 

1^1 


7 } sm a e sm a 

a + z — r —— , = a + 


1 S 1 

|S| 


1 — 277 cos a 


1 — e cos a 


§ 17. I here add another solution of Kepler’s problem in 
Series, the denominators of the terms being powers of 1 — e cos a. 

If 

x £= a + efx 

then, by Lagrange’s theorem, 

• ^^ a + ^’ a + Ma f2aWa + lh^ fSa Ffl + &e ‘ 


Now let 

therefore 

dx 1 


Far = 


dx 


I — efx da 


e 2 fafa e* d f 2 af"a e 4 d 2 f*af"a 0 

■ *■ - ~r; + - — jz - ' 7 T -+ & C . 


da 1—efa (1—efa ) 2 1.2 da (1 —efay 1.2.3 da 1 (1 — efay 

_ efa e 2 faf a ef d, f 2 af r a 

1 + 1 — efa + (1 — efay + 1.2 da (I — efa )' 1 + C * ’ 

whence, integrating, 

efa ( e 3 f 2 af"a e x d f 3 af tr a ( 

* ~~ ^ + I — e/'a + 1.2 (i - efay + 1.2.3 da (1 - efay *" C * 

In the case of Kepler’s problem 


fa =s sin a, fa - cos a, f r a = — sin a\ 


and therefore 


x — a + 


e sm a 




sm 3 a 


e* d sin 4 a 


1 — e cos a 1.2(1 —0 cos a ) 2 1.2.3 da (1— e con af 

e 5 e£ 2 sin 5 a 


1.2.3.4 cfo- (1 — g cos «) 


&c. 


= a + 


e sm a 
1 — e cos a 

sin 3 a 


2 (1 — e cos a ) 2 


-\2 


sin 3 a cos a 


e sin 5 a 


3 I (1 — e cos af 2 (i—e cos a ) 3 j ’ 
sin 3 a _ . sin 5 a 


..... w e sin 5 a cos a , e 2 sin 7 a ) 

-_ 2 C _22_*_- x 5 _( 

24(i — e cos «) 2 ^(i— e cos a ) 2 (1 — geos#) 3 (1—ecos«) 4 j J 


e 5 f 

-—1 20 


— &c. 

1877, August 4. 
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